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Abstract 



In this work the production of vector mesons in diffractive DIS is analyzed, with an em- 
phasis on the impact of the internal spin structure of a vector meson upon its virtual photo- 
production rate. For the first time, the full sets of s-channel helicity conserving and violating 
amplitudes were derived for pure S and D wave vector mesons. In the course of analytical 
and numerical investigation, we found a striking difference between S and D wave vector me- 
son production. First, in the case of helicity conserving amplitudes, D wave vector mesons 
were found to exhibit dramatically different Q^-behavior due to abnormally large higher twist 
contributions. Besides, we found that helicity conserving amplitudes for S and D states are 
numerically comparable, which makes the physical p meson very sensitive to D wave admix- 
ture. Second, the helicity violating amplitudes in light vector meson production turned out 
large both for S and D wave states. In the case of heavy quarkonia, these amplitudes became 
suppressed by the non-relativistic motion for S wave mesons but they were still large for D 
wave states. 

Aiming at producing sort of a manual on the diffractive vector meson production, we carried 
out all calculations in a pedagogical manner with maximally full presentation of intermediate 
steps and accompanied them with detailed qualitative discussion. 
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1 Introduction 



In the past 30 years the particle physics theory has proved numerous times to provide a 
good, consistent, unified description of the great variety of nuclear, low and high energy 
particle physics experiments. Being based on the ideas of QFT applicability, gauge approach 
to fundamental interactions, symmetry and naturalness considerations, the Standard Model 
managed to explain virtually all phenomena in electromagnetic, weak and strong interactions, 
to predict new particles and effects. Though the questions of fundamental origin lie beyond 
the scope of the Standard Model, its precision in description, for instance, QED phenomena 
reaches the magnitude of 10"^''. 

However, the current situation is not that optimistic in the domain of strong interactions. 
The gauge-based formulation — the Quantum Chromodynamics (QCD) — seems to offer 
reasonably good description only of the energetic enough processes (more accurately: only 
when every vertex involves at least one highly virtual particle) thanks to the famous asymptotic 
freedom. The major difficulty lies in the behavior of the QCD coupling constant: OsiQ"^) 
exhibits infrared growth and becomes comparable to unity at ~ 1 GeV^. The net result 
is that the perturbation theory — the only prolific universal treatment of various high-energy 
processes — fails to give even qualitative description of low-energy, essentially non-perturbative 
phenomena. Additional difficulties arise from the non-abelian nature of QCD, chiral symmetry 
breaking, non-trivial QCD vacuum, instantons etc. 

On the other hand, many separate concepts have been developed, which do not cling to 
the perturbative QCD (pQCD) and provide reasonably good description of phenomena in 
their applicability regions. The fundamental problem of the theory of strong interactions is 
that these heterogeneous approaches do not matchQ. They do not comprise a unified picture 
of strong interactions. Given such a lack of universal, rigorously derived results, one must 
admit that the subject of our investigation belongs to the realm of phenomenology rather 
than rigorous theory. 

1.1 DifFr active processes and Pomeron 

In the light of these problems, the careful examination of regions where two or more approaches 
overlap (or visa versa, where neither of the concepts exhausts the interaction) are of great 
interest. Diffractive Deep Inelastic Scattering (DIS) is exactly one of these fields. 

A typical hard DIS process (Fig.^) occurs when a virtual photon^] strikes a proton to 
produce a hard system X with large invariant massQ s and large enough multiplicity, final 
state hadrons being distributed over whole rapidity range approximately smoothly. 

However, as it was noted long ago, sometimes the proton survives, being only slightly 
defiected, and a virtual photon turns into a so-called diffractive system Xdifjr with invariant 
mass Mj^jrjj. <^ s. In this process the proton and the diffractive system are naturally separated 

^ Just a few examples of poor accordance among various approaches: the quark generated ladder diagrams 
do not appear to correspond uniquely to any of experimentally observed Regge trajectories. Another example 
is the vague status of as = const BFKL results in true QCD. 

^We will always imply that the virtual photon is emitted by an electron, which means the photon is always 
space-like: if q is photon momentum, then = —q^ > 0. 

^In hard DIS phenomenology this quantity is usually labeled as W^. However, for simplicity we will use 
notation s. 
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p p b) 

Figure 1: Examples of deep inelastic scattering process: (a) hard DIS and (h) diffractive DIS. 
In the latter case Mj-jj^, -C s and the process proceeds via pomeron t-channel exchange. 

by a large rapidity gap and a condition which appears necessary for the rapidity gap formation 
is -C s, or in terms of Bjorken x 

X = ^ « 1 . (1) 

s 



This is one of the most common cases of diffractive DIS (DDIS) processes. In fact, the 
class of diffractive processes is not confined within DIS; it is much broader. There are many 
other reactions which possess the generic features — the rapidity gap and smallness of 
— and therefore can be classified as diffractive processes (for a recent review see 

How can a typical diffractive process occur? Certainly, it must be kind of a peripheral 
interaction: if the photon struck directly one of the valence quarks, the proton would 'explode', 
providing no way for the large rapidity gap formation. What remains is the possibility of the t- 
channel exchange by not-too-energetic 'particle' (Fig.|l|b), which would be a natural mechanism 
of the experimentally observed weak proton deflection and small Mj-j-jr^.. Further experimental 
features suggest that this 'particle' should be chargeless and colorless, its interaction with other 
particles should be of strong (not EM or weak) nature, its 'propagation' should be independent 
of the specific process (7p, 77, pp, pp, etc), and it should be of spin 1 (due to approximately 
s-constant pp cross section). In the early 60s this 'particle' was dubbed Pomeron (symbol 
IP). 

Further properties come from combining the Regge picture and BFKL results with ex- 
perimental observations (for a detailed review of Regge theory see 0). They include, first 
of all, the asymptotic equality of total pp and pp cross sections (the Pomeranchuk theorem). 
Formulated long ago, it was experimentally verified only recently. Then, the Regge theory 
predicts the power-like s-dependence of the total pp cross section a oc s^Caip-i)^ which has also 
been experimentally observed, with intercept 6jp = ajp — 1 ^ 0.08. On the other hand, the 
BFKL equation 0, ^ succeeded in reproducing such power-like dependence in QCD, but in 
a simplified case as =const. In this approach the hard pomeron is treated as two reggeized 
gluons — an ansatz used currently in diffraction phenomenology with great success. However, 
the predictive power of the BFKL approach for the numerical value of the pomeron intercept 
is still limited and not all issues with sensitivity of the result to the infrared region have been 
understood. For further reading on pomerons, a topic very intriguing by itself, we refer to |^ . 



6 



1.2 Vector meson production in difFractive DIS 

There are several possible final states X in a typical diffractive DIS (DDIS) process 7*p —>■ Xp: 
system X can be a real photon, a qq continuum pair forming two jets or qq bound state, for 
example, a vector meson. Let us now focus specifically on exclusive vector meson production 
in diffractive DIS. This reaction has been studied extensively at fixed target DIS experiments 
at CERN and FNAL and more recently by the HI and ZEUS collaborations at HERA. 

In an off-forward scattering 7^^ the s-channel helicity flip amplitudes can be non- 

vanishing. Because of the well known quark helicity conservation in high energy QCD scatter- 
ing, such a helicity flip is possible only due to the internal motion and spin-angular momentum 
coupling of quarks in a vector meson. This issue was accurately analyzed only in very recent 
papers |^ , where it was shown that helicity non- conserving amplitudes are not negligible, as 
they had been thought before. Thus, as such, the helicity flip amplitudes would offer a great 
deal of unique information of internal constituent motion and spin-angular momentum struc- 
ture of vector mesons, unaccessible in other experiments. In addition, the vector meson decays 
are self-analyzing and the full set of helicity amplitudes can be measured experimentally. For 
unpolarized incident leptons, the angular distribution of decay products is parameterized in 
terms of 15 spin-density matrix elements, which can be calculated via 5 — 2 helicity conserving 
plus 3 helicity violating — basic helicity amplitudes |^ . 

Despite the great deal of theoretical work on vector meson production in diffractive DIS 

0, |lO|, the above issue of sensitivity to the spin-angular momentum coupling has not 
been addressed before. Namely, in a typical vector meson production calculation, a vector 
meson has been implicitly taken as 15* state and at the same time an unjustified ansatz was 
used for qq ^ V transition spinorial structure, namely, of u''yf^u ■ type. Being a mere 
analogy of qqj vertex, this ansatz in fact corresponds neither to pure S nor to pure D wave 
state but rather to their certain mixture. Only in [|n|] the cases of IS* and 25* vector mesons 
were compared and the necessity of similar calculation for D wave states was stressed. Such 
calculations however have been missing in literature until now. 

In addition to purely theoretical needs, there are more issues that call upon a thorough 
analysis of the D-wave effects. For instance, different spin properties of the S- and D-wave 
production may resolve the long standing problem of the Z)-wave vs. 25'-wave assignment for 
the p'(1480) and p'(1700) mesons (as well as the u' and 0' mesons). Furthermore, the deuteron 
which is a spin-1 ground state in the pn system is known to have a substantial D wave admix- 
ture, which mostly derives from the tensor forces induced by pion exchange between nucleons. 
Recently, there has been much discussion [|12| of the nonperturbative long-range pion exchange 
between light quarks and antiquarks in a vector meson, which is a natural source of the S-D 
mixing in the ground state p and u mesons. 

Motivated by these considerations, we undertook a derivation of a full set of helicity am- 
plitudes for diffractive electroproduction of pure S and D-wave qq systems for all flavors and 
small to moderate momentum transfer Aj^ within the diffraction cone. 

The structure of this thesis is following. In the next Chapter we develop a QCD-motivated 
formalism of pure S and D vector meson treatment. In Chapter 3 a thorough and detailed 
derivation of vector meson production amplitudes is given. In Chapter 4 we continue our 
analysis in the case of heavy quarkonia and discuss the most eye-catching differences between S 
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and D wave production amplitudes. In Chapter 5 we present some numerical results and finally 
we draw conclusions. In the course of calculations, we aimed at maximally full presentation of 
intermediate steps and results. Thus, we hope that the present work might also be helpful for 
anyone who intends to start his or her own calculation in this field. Some of the math-heavy 
calculations comprise the Appendix. 

We would like to underline that this work sheds a new light on the Deep Inelastic Scatter- 
ing as a means of hadron structure investigation. Indeed, one is accustomed to using DIS as 
a probe of the proton (i.e. target) spin structure. We demonstrate here that it is a powerful 
tool for studying internal spin structure of a vector meson as well. To this point, our work 
reflects only a few steps in this promising and still not fully explored domain. 



All the results given in this thesis were derived by the author. The only exception was the 
results shown in Fig.rU, which were calculated by I.Akushevich, Inst, of High Energy Phys., 
Misnk . The majority of the results derived in this thesis have been published in [|13[ and 
presented as talks in [111, [15 
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2 Description of a vector meson 



In this section we first introduce the vector meson fight cone wave function (LCWF) and show 
how it emerges in diagrammatic calculations. Then, describing S and D wave type vector 
particles, we give at once expressions for S and D wave vector meson spinorial structures, 
which we then prove by computing the normalization condition for LCWF. Finally, we also 
calculate V — >• e~^e~ decay constants to be used afterwards. 



While describing particle motion in non-relativistic Quantum Mechanics, one usually deals 
with a configuration space particle wave function, which is a good description because the 
number of particles is conserved. So, when one has a system of particles and shows that the 
wave function corresponding to their relative motion descreases at large relative distances at 
least exponentially, one can speak of a bound state. 

In Quantum Field Theory (QFT) this approach needs an update, since the field function 
becomes an operator in Fock space. Besides, since a bound state always implies the presence 
of interaction, the projection of a physical bound state onto the Fock space of free, non- 
interacting, plane-wave state vectors has a rather complicated structure: 



We emphasize that in this decomposition quarks and gluons are assumed free, i.e. on mass 
shell. Coefficients q can be called 'wave functions' of the given projection of a physical vector 
meson, with |cip being the probability of finding a vector meson in a given state. 



Let us now outline how a wave function of a bound state appears in the diagrammatic language. 

In the non-relativistic quantum mechanics, the two-particle bound state problem can be 
immediately reformulated as a problem for one particle of reduced mass /i, moving in the 
external potential. This reformulation allows one to split the wave function into two factors: 
the wave function of the motion of the composite particle as whole and the wave function 
corresponding to the internal motion of constituents. The former part factors out trivially, 
while the latter wave function obeys the following Schodinger equation 



Since the wave function ■0(r) and the interaction operator V{r) exhibit good infinity behavior, 
one can rewrite this equation in the momentum representation 



2.1 Bound states in QFT 



Vphys) = Co\qq) + ci\qqg) + C2\qqgg) + c^lqqqq) + ... 



(2) 



2.2 LCWF and vertex factor 




(3) 
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In this notation, this equation can be viewed as a homogeneous non-relativistic Bethe-Salpeter 
equation for the wave function ip{p) that describes the relative motion of constituents inside 
a composite particle. 
Let us now introduce 



T{p) ^(^-e) (5) 



Then Eq.(H) can be rewritten as 



m 

This equation has an absolutely straightforward diagrammatic interpretation (Fig.^. One sees 
that r(p) stands for bound-state — > constituents transition vertex, with p being the relative 
momentum of the constituents, factor 1 / [ — E] describes propagation of qq pair and V{k) 
stands for the interaction between constituents. Of course, the kinetic energy 7^ E, the 

total energy, which is in fact negative, so no pole arises in the propagator. 





Figure 2: The diagrammatic interpretation of the integral equation for vertex function r(p) 
(p is the relative constituents momentum). 

In the relativistic case, i.e. in QFT, it is not clear a priori whether the whole picture 
that involves somehow defined wave function and representation of the vector meson as free 
non-interacting constituents would work at all. So, in our approach we will be aiming at 
introducing an appropriately defined wave function and demonstrating that hard processes 
involving vector mesons can be expressed in terms of expectation values of qq amplitudes 
between wave functions, i.e. we intend to treat a hard process in a probabilistic, quantum 
mechanics-like manner. 

In the following we will show that this program succeeds. Namely, we will introduce the 
radial wave function of the qq state of a vector meson as 



(the angular dependence of the wave function will be treated separately, see Sect.^.4[) Here 
r(g) is the vertex factor, is the eigenvalue of the relativistic kinetic operator of the on 
mass shell qq state and rriy is eigenvalue of the total relativisitic Hamiltonian, which is of 
course equal to the mass of the vector meson squared. Then, during an accurate and honest 
analysis of a hard process Feynman diagrams, we will always make sure that wave function (|^) 
automatically appears in calculations and the rest looks the same as if both fermions were on 
mass shell. If we see that fermion virtualities modify the results, or if different Fock states mix 
during hard interactions of the vector meson, it would signal the invalidity of the free particle 
parton language and consequenly the breakdown of the whole approach. This restriction must 
always be taken into account when obtaining and interpreting the parton model-based results. 
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2.3 Light cone formalism 

The term "light cone approach" to high-energy process calculations can have different mean- 
ings. Some prefer to re-formulate the whole QFT within light cone dynamics, introduce 
light cone quantization and derive light cone Feynman rules (on Light Cone Field Theory see 



It, [13 )• However, one should keep in mind that even within the usual QFT the light cone 
formalism can be freely used as a means to greatly simplify intermediate calculations. This is 
exactly the way we will use it. 



It was noted long ago [|T^ that the calculation of a high energy collision is simplied if one 
decomposes all momenta in terms of light cone n+,n^ and transversal components, which we 
will always mark with the vector sign over a letter (so called Sudakov's decomposition): 

< = ^(1, 0, 1) ; n't = -^(1, 0, -1) ; (n+n_) = 1 , (n+n+) = (n_n_) = 

= p+n^ + p^n^ + P^ ', P^ = '^P+P- — P^ ■ (8) 

Indeed, imagine two high energy particles colliding with momenta p^ and respectively 
and equal masses m. Then one can choose such a frame of reference that in the Sudakov's 
decomposition 

= p+n^ + p-n^ + f?^ , 
= q+n^ + q-U^ + 

quantities p_ and g+ are large p_, g+ ^ m while p+ = {p^ + m?)/{2p_) and q- = {(f + 
m?) I (2g+) are small. The total energy squared of these two particles will be defined as 

s = 2g+p_ . (9) 

Note that our definition of s is somewhat different from the more familiar one (p + q)"^ by terms 
oc m?,fp'. However, it is not of any imporatnce for us, since in the course of all calculations 
we will keep track only of the highest power s contributions, i.e. everything will be calculated 
in the leading power s approximation. 

ki 



k2^ 

Figure 3: Kinematics of V —>■ qq vertex on the light cone. Vector meson momentum qy is 
taken incoming, constituents momenta are outgoing. 

Let us now go further and examine the kinematics of a typical qqV vertex (Fig.^. The 
Sudakov's decomposition of all momenta reads: 

q'^ = qv+n'l + qv-n_ ; 

ki = ki^n'^ + ki^nt + = zqy+n^ + yqv-nt + ; 

k^ = A;2+< + k2-n>t - fc^ = (1 - z)qv+n>'+ + (1 - y)qv-n''_ - P (10) 
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so that 

ql = 2qv+qv- = ; k'^ + k!^ = q!^; k} ^ , (11) 

i.e. quarks can be off mass shell. Now let us introduce momenta k*, which would correspond 
to on mass shell fermions. The only component in ki subject to modification is component, 
or absolutely equivalently, the energy. Large components are insensitive to (reasonable) 
quark virtuality variations. So, to obtain the on mass shell momenta, one has to replace 

i._.^l±i!^t._ = !!!!±^. (12) 



Then the 4-vector 
squared is equal to 



= kl'' + kl^ (13) 



= = 2g+g_ = 2qy^ {kl_ + kl_) = . (14) 

And again we emphasize that the Feynman invariant mass (i.e. the total 4-momentum 
squared) of the virtual quark- antiquark pair is my. The quantity is the invariant mass of 
the free, non-interacting qq state (see |^). However, it is precisely M, not my that will govern 
the hard interaction of qq pair with gluons. 

Finally, it is useful to introduce the relative momentum of free qq system: 

2p, = {ki-k;),. (15) 

Then, trivial algebra leads to 

M2 = 4m2 + 4p2; / = -p2; {pq) = 0, (16) 

where p is the 3-dimensional relative momentum in the qq pair rest frame of reference. Its 
components are 

P = {p,Pz); P = k; p, = ^{2z-l)M. (17) 
2.4 Spin structure of a vector particle 

Let us start with a well known example of a deuteron, which is a non-relativistic analogy of a 
vector meson: they are both vector particles built up of two fermions. To have the correct P- 
parity, proton and neutron must sit in the spin-triplet state, thus leaving us with two possible 
values of their angular momenta: L = and 2. 

In the conventional non-relativistic language one describes the spin-angular coupling by 
the Clebsh-Gordan technique. The non-relativistic Feynman diagram calculations can be best 
performed in an alternative approach. Here a deuteron, being a vector particle, is described by 
a 3 dimensional polarization vector V. So, while calculating high energy processes involving 
d pn transitions, one can use the following spin structure of deuteron-nucleon-nucleon 
vertex: 

0;^r0,-v (18) 
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Since both nucleons can be treated on mass shell, only two terms enter Fj, which can be 
written as: 

0+ [u{p)a' + w{p){3py - ^V)^-"] <Pp ■ V' (19) 

Here a* are Pauli matrices and p is the relative proton-neutron momentum. One immediately 
recognizes here spin structures corresponding to pn pair sitting in S and D waves respectively. 
In particular, squaring the above expression gives 

(VV*) for 

3(pV)(pV*) - (VV*)p^ for SD interference 

3p2(pV)(pV*) + (VV*)p^ for (20) 

Now, let us go relativistic and turn to vector mesons. The polarization state of a vector 
particle is described by a four- vector V^. Therefore, a general form of qqV vertex has the form 

uT^u ■Vf,-T{p) , 

where r(p) is the familiar vertex factor. Up to now, it has been customary in literature to 
choose the simplest form of the spinorial structure F^: 

u'^^u ■ ■ T{p) . (21) 

However, one must admit that ( |2l| ) is simply an analogy of gg7 vertex and does not reflect the 



true internal structure of a vector meson. It is known that the correct spinorial structure 
corresponding to pure the 5* wave qq state reads 

^^ = ^^- mT^ = - ^(m'+L)) " '^-^^ • ^22) 

It is implied here that spinorial structures are inserted between on mass shell spinors in 
accordance with our principal guideline (see discussion in Sect. p.2|) . 

Once S wave spinorial structure is established, the expression for D wave can be obtained 
by contracting 5* wave with the symmetric traceless tensor of rank two SpiPj — Sijp"^, rewritten 
in the Lorenz notation. To do so, one should replace 

■ A - M Mil 

Pi ~^ Pfi 1 '^ij ~^ 9tiv "T 

(in the qq pair rest frame of reference q^ = (M, 0, 0, 0)). However, since q^ inserted between 
on mass shell spinors gives zero due to the Ward identity, one obtains the required tensor in 
the form 3p^p,^ + g^uP"^- Its contraction with yields 

2\r- 2 / \ ( 1 M + m \ 
= {3p^Pu+gf,uP )-<Sup7p = P 7^+(M+m)p^ = ( p g,,^ H — — p^p^ \ Iv = V^y^i^ . (23) 



We will prove below that structures (|2^), ( pBf ) after being squared indeed perfectly reproduce 
(pOD, i.e. they indeed correspond to pure S and D waves. 



The quantities S^^ and V^^ used in (p2D, (^3]) have the meaning of S/D wave projectors, 
which will be used in all subsequent calculations. Namely, all calculations will be at first 
performed for the naive qqV vertex (^) and then we will apply the projector technique to 
obtain expressions for S and D wave states. 
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2.5 Vector meson LCWF normalization 



Before tackling the diffractive vector meson production process, we first should have a pre- 
scription of normalization of the vector meson wave function. 



2.5.1 Naive qqV vertex 




Figure 4: Diagram used for normalizing the vector meson LCWF. The amplitude of this 
diagrams is set equal to 2q^i. 

A natural way to normalize the wave function of a composite system is to put the ampli- 
tude given by this diagram in Fig.^ equal to 2g+i. Here extra leg carries zero momentum but 
couples to the fermion line as 7^n^. Note that for a charged composite particle (a deuteron) 
this is precisely setting the electric for mf actor equal to unity in the soft photon limit. 

As described above, we first treat qqV vertex as u''~^^u ■ In this case the general 

expression for this amplitude is 

(—1) r 4 Sp{iViT ■ i{k — gy + m) ■ iV^T* ■ i{k + m) ■ ih_ ■ i{k + m)} 



(27r)^ J [A;^ — + ie] ■ [/c^ — + ie] ■ [{k — qvY ~ '^'^ + ^c] 



{271)^ J [k"^ - m? + itf[{k - qvY - m? + ie] 

where A^c = 3 is a trivial color factor originating from the quark loop. We deliberately 
recognized Vi and V2 as distinct entities just to make sure later that such a loop is indeed 
diagonal in polarization states. 

The first step is to rewrite this expression in terms of Sudakov's variables. As usual, one 
implements decomposition 



k^" = zqv+n1 + yqv~n^- + k^ ; qy = 2qv+qv 



m 



V 



and transforms 



d^k = -rriycfkdydz . 



Now we note that vertex functions F do not depend on y (and neither does the trace, as will 
be shown later), so we can immediately perform the integrations over y by means of Cauchy 
theorem. Indeed, since the integral 

DO 1 1 

dy ^ ^ , (25) 

00 [yzTUy — (/c2 -|- 77^2^ _|_ ^gj2 ^(^i _ — z)my — (k"^ + m?) + ie] 
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is convergent and has good infinity behavior, one can close the integration contour in the most 
convenient way. To do so, one should analyze the position of all poles on the complex y plane: 



k"^ + ie k"^ + m? ie 
^1.2 = ——2 — ^; 2/3 = 1- + 



ziTLy zruy' {1 — z)my {1 — z)my 



One sees that if ^ < or 2 > 1, a// poles lie on the same side of the real axis in the complex 
y plane, which leads to zero contribution. The contribution that survives comes from region 
< 2; < 1, which has in fact a simple physical meaning: all constituents must move in the 
same direction. In this region, we close the integration contour through the upper half-plane 
and take a residue at y = 2/3. Physically, it corresponds to putting the antiquark on mass 
shell. After this procedure, one gets for (^): 



27ri 



2m (1-z 



i2 



1 - z)ml [yzm\r - (F + m2) + ie]^ y=yi (1 - z)m\ [p + m? - z{l - z)ml.f 

2m 1 1 
~mlr z'^{l-z) [M^-ml] 
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One immediately recognizes here the same two-particle propagator as in (|^). Therefore, the 
equation for the amplitude reads 

A = ^J^, ■ I d'k-J^ ■ IV^P • (-lsp{...}] . (26) 



(27r)3 J z^{l-z) V 2 

Now we turn to the trace calculation. To do this, we will first show that thanks to the 
presence of n_ vertex, both constituents indeed can be treated in the numerator algebra as if 
they were on mass shell. Indeed, applying the Sudakov decomposition to the 7 matrix, one 
has 

7^ = j+n1 + -f.n't + 7^ ; 

7+ = ^- = ^(^0 + ^3) , 7- = = ^(70 - 73) • 

Now let's decompose the propagator numerator of the constituent, to which this n_ leg couples: 

k + m = A;+7_ + A;_7+ — k^ + m (27) 



and rewrite it using notation of (|T0D,(|T2D as 

k + m = A;+7_ + /c* 7+ — k^ + m + (A;_ — fc* )7+ = k* + m -\ — 7+ . (28) 

In other words, we expressed the virtual quark propagator as the sum of on shell quark 
propagator and an additional "instantaneous interaction" term. However, since n_ is inserted 
between two {k + m) factors, this item does not work due to identity 7+7+ = 0. The net 
result is that wherever h_ appears, both constituents can be treated on mass shell in the trace 
calculation. This observation is very important for us, since, as we remember, it validates the 
whole approach. 
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Having established this property, we can now easily calculate the trace. In our case, the 
easiest way is to do it covariantly, without involving further the Sudakov technique. Since 
quarks in numerator can be treated on mass shell, we first note that 

[k + m)n-{k + m) = 2{k*nJ){k* + m) 

so that 

-\sp{...} = -zq+Sp{V,ik* - q + m)V2*Ck* + m)} = -2^g+ [m\V,V*) + A{V^p){V*p) 



where p is the relative quark-antiquark momentum [see (16)]. Note that in the antiquark 
propagator we replaced k — qv —>■ k* — q, since the antiquark is now put on mass shell. 
Besides, we explicitly used here gauge condition (qV) = 0, which means that polarization 
vectors must be written for on mass shell qq pair, not the vector meson, — another important 
consequence of our approach. Substituting this into (p6D, one gets 



(2vr) 



/ "^'^IfT^''^'' [-M\V,V*) - A{V,p){V*p)] . (29) 



A prominent feature of this equation is the orthogonality of Vl and Vr polarization states — 
the necessary condition for any normalization prescription. 

The next step is to realize that the integral can be cast in the form of c/^p integration by 
means of ,^4 4 

-(fk = —dpzd^p = -TJ'^^P ■ 



z{l-z) M M 

Tus, the final expression for normalization condition is 

N r A 

1 = y d^^-m-M\v,v;) - A{y,p){y;p)] . (so) 

We see that the expression being integrated is explicitly spherically non-symmetric, which is 
a manifestation of a certain D wave admixture. Thus we now apply projector technique to 
obtain results for S and D wave states. 

2.5.2 Normalization for S wave vector meson 

The correct expressions for pure S/D type vertices can be readily obtained with the aid of 
projector technique. Namely, to obtain an expression for S wave, replace 

. (31) 

Such a replacement for Vi leads to 

-M\V,V;) - A(y,p){V;p) ^ -M\V,V;) - ^^M_iy^p)(y;p) 
Then, one applies the same replacement to V2 to obtain 
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Therefore, the answer for S wave states reads 




(32) 



which is manifestly spherically symmetric. 



2.5.3 Normalization for D wave vector meson 

Results for D wave states are derived in the same way. The replacements 

- M^p\ViV;) + 3M^p\Vip){V*p) . 

After angular averaging 



V^V^^ lead to 



(33) 



1 



(PiPj) ^ 3P 

one gets the normalization formula for D wave state: 




(34) 

Several remarks are in order. First, S wave D wave transitions are forbidden. Indeed, 
such an amplitude will be proportional to 



(35) 



which vanishes after angular integration. Then, we emphasize that the structure of results 
(|32D , (0), (|35|) is absolutely identical to Eq.(|^). This fact can be viewed as the proof t\iaX 
spinorial structures ([2^), (PB|) indeed correspond to pure S and D wave states. 



2.6 Decay constant 




1/ 



q-k 



0-) 




6) 



Figure 5: Normalizing LCWF to r{V — e+e ) decay width: (a) the diagram for y — 7* 
transition, (b) the diagram for V e^e~ decay. 

An additional normalization condition consists in relating the vector meson wave function 
to the experimentally measurable physical quantity — V e^e~ decay width (Fig.|^). The 
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loop at Fig.^ describes transition V ^ Y and enters the amplitude of the decay V e+e 
(Fig.pj). Let us define the decay constant via relation 



(36) 



So defined fv has dimension dim[/y] = m^. The quantity cy refiects the fiavor content of a 
vector meson (in the previous calculations it simply gave unity) and is equal to 



1 



1 2 



cv 



v^'3v^' 3' 3 



(37) 



for p, uj, (j), J /if} mesons correspondingly. 

Knowing that such a loop does not mix polarization states, we can multiply both sides of 
eq.(^) by V* and get the expression 



V 



(2vr) 



-Nr ■ / (fk 



Sp{iV* ■ i{k + m) ■ iVT ■ i{k — qv + m)} 

(A;2 — m? + ie) ■ {{k — qyY — fn"^ + ^e) 



(38) 



Calculations similar to the above normalization condition derivation yield (for the naive 
type of vertex) 



/v 



Nr. 



(2vr) 



dz 



zil-z) 



(fk ^iJv[-M'^iyv*) - A{yp){y*p)] . 



(39) 



Applying now the projector technique, one gets in the case of 5* states (after proper angular 
averaging) 



(40) 




and in the case of D wave states 



/ 



Nr. 



J P -ipD- 



32 



3 M + 2m 



Finally, one can write down the expression for the decay width in terms of fv'- 



e e 



4:Tr\A\ 



3my 



2 2 



(41) 



(42) 



This formula can be used to extract the numerical value of fv from experimental data. 



2.7 Ansatz for LCWF 

In order to perform some model evaluations, one should fix a particular form for the vector 
meson wave function. We will give just two examples often used in calculations: Coulomb-type 
and oscillator type wave functions. 
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2.7.1 Coulomb wave functions 



As one can see from Eqs.(|32D and (@), a LCWF should contain M dependence in form of 
M~^/^ (the normalization condition for non-relativistic wave function deals with / cPp\'4'nr\^ )• 
So, in this ansatz the normalized wave functions read 



Cl 




1 












-al) 






+ alf 


cd 




1 




(p2 





Mp') = ^7::^^7y^- (43) 



with normalization constants to be determined from Eqs.(P^ and (0). Here parameters 
are connected to the size of a bound system: in the coordinate representation ipi oc exp(— rctj). 
For usual Coulomb functions Sa^ = = ai = 1/RBohr, where Rsohr is the Bohr radius. 
However, this relation should be treated with care in the case of qq quarkonia, where the 
quark-antiquark potential is quite complicated and therefore ctj should rather be considered 
as free parameters. 

As an example, let's compute the ratio fo/fis with these wave functions in the non- 
relativistic limit, when M ^ my ~ 2m. It will be useful later in the analysis of heavy vector 
mesons. 

fD 




fs ^ 9 \ Hh mlRl ■ ^^^^ 



I'D 

Here i? = l/a is a typical size of a vector meson in a given state. 

2.7.2 Oscillator type LCWF 

In this ansatz one has 



'^2 A 2__2d2\ / P^-^i 



^2. = ^1 - -p i?, j exp 

'^D _ ( P^^D 



Note again that for purely oscillator potential one also has relation Ra = Ri = R2, which 
might not hold in our case, since the oscillator type potential is also a crude approximation of 
the true quark-antiquark interaction. A similar calculation for fo / fs ratio yields 



fn 160 



Rl 1 

^ (46) 



fs ^ 3 \ Rl mlRl' 

Note that the oscillator ansatz gives at least times larger fo/fs ratio. Therefore, one must 
admit that such model dependent calculations are plagued by a great deal of arbitrariness, 
which will inevitably affect the results. 
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3 Vector meson production amplitudes 
3.1 Preliminary notes 




p p p p 

a) b) 



Figure 6: The QCD-inspired diagrams for 7*p Vp process with two gluon t-channeL Only 
Diagr. (a) does contribute to the imaginary part of amphtudes. 

Having set up the notation and defined and described a vector meson by itself, we are now 
ready to evaluate the full set of amplitudes of its off-forward virtual diffractive photoproduc- 
tion. 

In the pQCD motivated approach to this process the pomeron exchange is viewed as a 
two-gluon exchange as it is shown in Fig.^. Using the scalarization procedure, we will split 
the diagram into 2 pieces and will treat each of them separately. The upper blob describes 
the pomeron-assisted transition of the virtual photon into a vector meson. In the perturbative 

2 

QCD approach, which is legitimate here due to the presence of the relevant hard scale Q = 
+ z{l — z)Q'^, the qq fluctuation of the virtual photon interacts with two hard gluons and 
then fuses to produce a vector meson. This interaction is described by four diagrams given in 
Fig.0, with all possible two-gluons attachments to qq pair taken into account. All of them are 
equally important and needed for maintaining gauge-invariance and color transparency. The 
latter property means that in the case of very soft gluons the upper blob must yield zero, for 
the qq pair is colorless. 

The lower part of the general diagram Fig.^ is of course not computable in pQCD. The 
physically meaningful procedure is to relate it to the experimentally measurable gluon density. 
To do so, we will first calculate this lower blob in the Born approximation and then give a 
prescription how to introduce the unintergrated gluon density. In the course of this procedure, 
we will neglect in the intermediate calculations proton off-forwardness and take it into account 
only at the very end, as a certain factor to the unintegrated gluon density. 

3.2 Notation and helicity amplitudes 

In our calculation we will use the following Sudakov's decomposition (see also Fig.||a) 

kf, = -yp^I + zq^' + ; 

= + crq^ + (47) 

Here k and k are momenta that circulate in the quark and gluon loops respectively; A is the 
momentum transfer. Vectors p^' and g^' denote the light-cone momenta: they are related to 
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Figure 7: The content of the upper blob in Fig.^ in the pQCD approach. The true vector 
meson internal structure is approximated by qq Fock state. 




AAAA 




P-A 



a) 



b) 



Figure 8: (a): A particularly useful convention of the loop momenta (only transverse com- 
ponents of quark momenta are shown), (b): the propagator notation used while calculating 
denominators. The crosses denote on mass shell particles. 



the proton and virtual photon momenta as 

P^=Pm' + ^?/; gM = ?/-^V; g"=p" = 0; x = ^^l; s = 2{p'q'). (48) 

As it was mentioned in the Introduction, the condition x -C 1 is necessary to speak about 
diffractive processes. The longitudinal momentum transfer can be readily found from kine- 
matics (see Fig.||). To the higher power s terms it reads 



ml = p^ = {p- A)2 =ml- 2{pA) + A^; a 



s 



ml = {q + Ay = -Q' + 2{qA) + A' ^ = x + ""'^ + . (49) 

The final vector meson momentum reads: 

qv, = g; + ""'^^^' V + A, . (50) 

Finally, throughout the text transverse momenta will be marked by the vector sign as k and 
3D vectors will be written in bold. 
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There are several possible helicity amplitudes in the transition 7^^ Vxv- First of all, 
both photon and vector meson can be transversely polarized. The polarization vectors are 

eT, = e^; ^TM = ^M + ^^(p'-g')M- (51) 

Note that we took into account the fact that the vector meson momentum has finite transverse 
component A. Then, the virtual photon can have the scalar polarization (which is often called 
longitudinal; we will use both terms) with polarization vector 

eo^, = ^{q +xp')f,. (52) 

Finally, the longitudinal polarization state of a vector meson is described by 

A2-M2 




W + ^f^j ■ (53) 

Note that, as we already mentioned, in the self-consistent approach we must take the running 
polarization vector for the longitudinal polarization state. It depends on M, not my, which 
reflects the fact that in our approach we first calculate the production of an on-shell qq pair 
with (whose dynamics is governed by M) and then projects it onto the physical vector meson. 
We stress that this projection will automatically arise in the course of usual Feynman diagram 
evaluation. 

Thus, there are 5 different amplitudes: 



Ay) 



-Ay) (54) 

The first two are helicity conserving amplitudes. They are dominant and almost insensitive 
to the momentum transfer A. The next two are single helicity flipping amplitudes. They are 
unavoidably proportional to | A| in the combination (eA) or {V*A) and would be vanishing for 
the strictly forward scattering. Finally, the last amplitude corresponds to the double helicity 
flip and will be proportional to {eA){V*A). 
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-^T (A 



3.3 General amplitude 

We will take diagr.(c) at Fig.|^ (it is shown in fig.§) as a generic diagram and perform a 
thorough analysis for it. It turns out that the other diagrams are calculated in the same 
fashion. 

The general expression for the amplitude given by diagr.(c) reads: 



iA 



{2tiY J (271) 



p - kl+TUp 



(p 



+ ie 



Hgi^'t^')u, 
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fcf — /i^ + le ^2 — /i + It 



Sp |iee i{k4 + m) {—ig^y'^t^) i{k^ + m)iV* i{k2 + ^) {—igi^t^) i{ki + m)| 
— + ie] [/cl — ""^^ + ^e] [k'i — m? + ie] [k1 — m? + ie] 

(55) 



Here cy is the same as in (pTf) and F is the famihar qq ^ V vertex function. Note that we intro- 
duced 'gluon mass' fi in gluon propagators to account for confinement at a phenomenological 
level. 

Let's first calculate the numerator. 
3.4 Color factor 

If we consider strictly forward gluon scattering off a single quark, we have 

1 1111 N'^ — 1 1 9 

j^MtV] . S,,.S,,,Sp{tV} = j^y.s-^S,s = 5^ = jC. = - (66) 

However, we should take into account that quarks are sitting inside a colorless proton, whose 
color structure is 

i^coior = -^e-'^q'^q'q' (57) 

In this case there are two ways a pair of gluons can couple 3 quark lines (see Fig.|^). In the 
first way both gluons couple to the same quark. Since the quark momentum does not change 
after these two interactions, the nucleon stays in the same state: (A^|A^) = 1. In the second 
case gluon legs are attached to different quark lines, so that extra momentum k circulates 
between quarks, which gives rise to the factor (A^| exp(mri — iKr2)\N), i.e. to the two-body 
formfactor. Therefore, for the lower line instead of 

one has 



a'b'c' 



^e"'"' {36aa'Swti„tX' + 6Saa'ttb'tZ'{N\ exp(^/cri - inr2)\N)) e 
= Sp{t^t^} - Sp{t^t^}{N\ exp{iKri - iKr2)\N) 

= ^6ABil-{N\exp{iKr^-iKr2)\N)). (59) 

Note also that a similar calculation for Nc number of colors would yield the same result. Thus, 
the overall color factor is 

^CfNcV{k) = 2V{k) = ^CfN,{1 - {N\ exp(mri - iKr2)\N)). (60) 
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Figure 9: The ways two glons can couple a colorless nucleon. 



3.5 Scalarization of upper and lower parts 

As known, the highest power s contribution comes from so-called nonsense components of 
gluon propagator (density matrix) decomposition: 



The lower (proton) line gives then 

u'{p) ■ q'{p — ki + mp)q' ■ u{p — A) 
To the highest power s order, it can be rewritten as 



u'^q'p'q'up 



= 7;Sp{p'q'p'q'} = . 

forward ^ 



(61) 



(62) 



(63) 



As we mentioned, the effect of off-forwardness (skewedness) will be taken into account later. 
So, combining all factors, one has for numerator of Eq.(|55|) 



(47ra,)2 V^TTaem cv ■ -CfN^V{k)—s'^ ■ Sp [e {h + m) q' {h + m) V* (/cs + m) q' {ki + m)| 



= (47r«,)2 y/4^i^ cv ■ 2CFN,ViK) ■ 2s^ ■ V) . 

Note that we factored out 2s^ from the trace because it will appear later in all trace calcula- 
tions. So, the resulting expression for amplitude (0) looks like 



(64) 



(27r)4y (27r)4 [(^ _ /,^)2 _ ^2 + ,g 

r*/W(7* ^ V) 



[ki — im? + ie] [fc^ — m? + ie] [ki — m? + ie] [ki — rn^ + ie] 

One can now immediately write similar expressions for the other three diagrams (Fig.|^ a,b,d). 
Indeed, they will differ from Eq. ( |65D only by the last line: they will have different expressions 
for traces and propagator structures (i.e. the exact values for /cj.). 



(65) 



3.6 Denominator evaluation 

Now we turn to the calculation of denominators. As usual, we implement Sudakov's decom- 
position (^) and make use of relation 

d^k = -sdy dz (Pk . 
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The physical picture of the way we will do the resulting integrals is the following. We are 
interested only in the imaginary part of these diagrams. In fact, it can be shown that at the 
level of accuracy used here the diagram in Fig.^a gives rise only to the imaginary part of the 
amplitude. The real part is given by Fig.^j and can be readily found from analiticity (so that 
here is no need for additional calculations), however it turns out small due to smallness of 
pomeron intercept, so we will neglect it in our subsequent calculations. 

The imaginary part is computed by setting three particles in the s-channel cut on mass 
shell (this is illustrated in Fig.||b). One way to do so is to apply Cutkosky rule to modify 
our expression. Another, more straightforward and 'honest' way is to calculate three of the 
integrals (namely, over a,l3,y) via residues. That's what we are going to do. 

The details of this calculations are given in Appendix 0. Here we cite the result: 

r 

dy dz da d(3 



[all propagators] 

2Tri\^ f dz ^ . 1 1 / N 

-) ■ J 1(1— y'Av('-. * ) • |j:j + ,„. + ,(i_,)Q.,(,,2 + ^.,. (66) 

Here ki is the transverse momentum flowing through photon vertex along the fermion line. 
Particularly, for diagr.(c) it is equal to ki = k — {1 — z)A — k,2 (with the specific quark loop 
momentum choice given at Fig.(|^)). 

Let us now stop for a moment and take a look at this result. As we know, aside from longi- 
tudinal momenta integration, we also have outer integrations over transverse gluon momenta 
/g°° d'^K,. It will also be shown later that the upper loop will always (except for double spin flip 
amplitude) give factor k'^ in numerator. Thus, we will deal with integral 



d K-r^ ^ . (67) 

{K^ + fj.''y[(ji-k + {l-z)Ay + m^ + z{l-z)Q^] 

This integral is naturally split up into three regions: < jj^, > > n"^, and > . 
(Here we marked all low energy scales, namely, /i^, m^, by a general parameter /i^.) A crucial 
observation here is that the largest contribution comes from the perturbative intermediate 
region 

< < g2 ^gg^ 

due to an extra large logarithmic factor \og{Q / ji^) (which of course implies Q ^ jj?). It is 
precisely the presence of this large logarithm that justifies our perturbative QCD 
approach. Thus, we now established that in our problem the pQCD approach is valid only 
in the leading log approximation and has only logarithmic accuracy. 

Thus, the amplitude for diagr.(c) has the form 



n 2 ^ ^ ( 2vri \ =^ 1 

A = V^T^Oiem 4:CfNcS Cy ' -S -S ■ [ • — — 

Z Z \ S / [ZlT) 



-^^d^Uv{z,l?) [ 7^^^^ (47ra,)2 — . (69) 

zil-z) ^"^^ ' 'J (K^ + fi^y"^ [kl + m'' + z{l-z)Q^] ^ ' 
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After bringing all coefficients together, one gets 

dz 



IS- 



CpNcCy^JA-Kae 



(27r)^ 



zil - z 



(PkV{k) 

{k^ + /i2)2 



• (70) 

A;f + m2 + z{l - z)Q^ 

The other diagrams are calculated in the same way. The most important difference is that 
for each diagram we will have a propagator l/[kl + -\- z{l — z)Q'^] with its own definition 
of ki, the transverse momentum in photon vertex: 

diagr.a kia = k — {1 — z)A = f— - A 

diagr.b kn, = k — {1 — z)A + k + -A = f+ k 



diagr.c k 



Ic 



'1 — z)/S. — — A = f— K 
' 2 



diagr.d k\^ = A; + zA = f+ -A 



(71) 



Here r = A; — (1 — 2^;) A/2. Thus, the whole expression for the imaginary part of the amplitude 
is 



A 



is- 



(27r)^ 



dz 



zil - z 



-d'^k'ipviz, k'^ 



d'^KV{K) 
?2 I n2\2 



l-z 



z kl^ + m? + z{l- z)Q'^ kff, + + z{l - z)Q^ 



+ 



kl + m2 + zil - z)Q^ l-zkl^ + m? + z{l - z)Q^ 
3.7 Final exact results for the naive vertex 



(72) 



The only thing left to be computed is integrands /*(7 V). For convenience, their calculation 
is also given in Appendix p. It turns out that the results can be written in the same way for 
all four diagrams via ki given by (O); i-^- all quantities: 

l-z 



IjW j{h) j(c) 



l-z 



.jid) 



can be written in a similar way: 



T 



{eV*){m^ + kki) + {V*k){eki){l - 2zf - {ek^Vk^) 
L^L -Az^{l - zfQM 
T^L 2z{l- z)M{eki){l-2z) 
L^T -2z{l - z)Q{l - 2z){V*k) 

Therefore, we can cast amplitude (^) in a compact form with aid of functions $i and $2: 
111 1 



$9 



[r + KY + Q' {r-Kf + Q' (f + A/2)2 + (f-A/2)2 + g' 



(73) 
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r + K f-K f+A/2 r-A/2 , 

•^1 = =2 =2 + — =2 + — =2 74 

{r + Rf + O' (f-/?)2 + Q' (f+A/2)2 + g' (f-A/2)2 + g' 



and 



With these functions, for the naive qqV vertex the whole expression in square brackets in (ff^ ) 
wit sign minus (which we denote here as I\^Xy) has the form: 

III = -AQMz\l-zf^2] 

Itt = (ey*)[m'$2 + (^$i)] + (l-22)2(fcV*)(e$i)-(a)(r*$i); 
Itl = 2Mz(l-2)(l-2z)(e$i); 

Ilt = -2Qz{l-z){l-2z)iVk)^2. (75) 
3.8 Gluon density 

We have just calculated our reaction in the Born approximation. We understand however 
that the nature of t-channel exchange in diffractive processes is much more complicated. A 
physically meaningful procedure is to related this and other diffractive processes to the ex- 
perimentally measurable gluon density. Here we explain a prescription of how to introduce 
unintegrated gluon density in answers calculated in the Born approximation. 

A well known QED formula for Weizsaker- Williams approximations unambiguously defines 
the number of photons in an electron (in the Born approximation and to the logarithmic 
accuracy) 

an, - 2^ [i + U /5J J ^ ~ ^ -,2 ^ . UDJ 
Here (3 is the fraction of electron energy carried by photon. Quantity 

= V = ^^^^ 

is called photon /3-density. Thus, the unintegrated photon density is given by 

^(Born) ^ _^P2_ = ^ . (78) 
a log 77 

In QCD, in the Born approximation, one can define a similar quantity, namely, the uninte- 
grated gluon density in a proton, which after careful accounting for all color factors reads: 

= -^^^ = CpK^V{k) . (79) 

O log 77 

Therefore, a prescription of how to include unintegrated gluon density is as follows: 

77 

Finally, one has now to account for skewedness of this gluon density. Within diffraction cone 
(i.e. for A^ up to several GeV^) it is usually parameterized as 

Hx, ^, A) = exp(-ii?3iP A^) . (81) 

o log K"^ 2 

where B^jp ~ 6 GeV~^ is the so-called diffraction cone slope. 

Therefore, a general amplitude of transition 7^ — Vx^ has the following form 
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A{x,Q,A) = -is 



27r2 



Jo z{l — z) J 



I 



{l-z) 



V), 



(82) 



3.9 Final results for 5* and D wave amplitudes 

Now we can use projector technique to obtain results for S/D wave states. 

2Pa^P- _ rs_r , 2(Vp) 



m{M + 2m) 



m(M + 2m) 



(83) 



V- 



(M + m)(Vp) 



m m 
Note that (Vp) is 3D scalar product. While contracting, we encounter terms proportional to 
P/i ® 7/i which should be understood as 



7/. = WAV ^ P} + WA^ = V,^p, = -(2z - 1)M}. 



The result of this substitution reads: 
for ex' 



1 



It^t{V ^p} + It^l{1 = V,^p, = -{2z- 1)M} 



= m 



{ek)^2 - {e^i){l - 2zf 



for cq: 



IlMV ^P} + h^L{l = V,^p, = -{2z- 1)M} 



= -2Qz{l - z)i2z-l)m^^2 
So, the resulting integrands for S wave type mesons are 



(84) 



(85) 



(86) 



^L-*L 



-AQMz'^{l - zf 



1 + 



(1 - 2zf 2m 
4^(1 - ^) M + 2m 



(eV*)[m'$2 + (^^i)] + (1 - 2zf{W*){e^^ 



M 



M + 2m 



2m 



M + 2m 



{ke){kV*)^2; 



2Mz{l-z){l-2z){e<^?i 



1 + 



'1 - 2z) 



'.III 



4:z{l-z)M + 2m 



Mm 
M + 2m 



-2Qz{l- z){l-2z){V*k) 



M 



M + 2m 



'l-2z)iek)^2; 



(87) 
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and for D wave type mesons are 



(ey*)p^K<l>2 + (A;$i)] + (1 - 2zy{p' + + Mm){kV*){e^ 
-p2(a)(r*$i) - m(M + m){ke}{kV*)^2 ; 



-M(l - 2z) 



(e$i) ( k'^ - —pI ) + m(M + m){ek)^2 



-2Qz{l - z){l- 2z){V*k){Yy^ + + Mm)$2 . 



Equations (i^); (H) together with expression (|82|) constitute the uhimate sets of all helicity 
amplitudes. They give explicit answers for the vector meson production amplitudes within 
leading-log-approximation. 
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4 Analysis for heavy quarkonia 



The general answers (P7|), ( |55D are of course incomprehensible at a quick glance. Therefore, 
a further analysis is needed to grasp the most vivid features of the results and to disentangle 
s-channel helicity conserving and double helicity flip amplitudes. 

Since in the heavy vector mesons quarks can be treated non-relativistically, further sim- 
plifications in analytical formulas (^), (|^) are possible due to the presence of an additional 
small parameter p^/m^. 

In what follows we will first perform the twist expansion and then relate simplified ampli- 
tudes to the decay constants (^OD, (|iT|). We will then analyze twist hierarchy of the amplitudes 
and compare results for S vs. D wave states. Though we perform this analysis for heavy 
mesons, we wish to stress that all qualitative features {S vs. D difference, dependence 
etc.) will hold for light quarkonia as well. 



4.1 Twist expansion 

Here we are going to expand the amplitudes (or to be more exact, the quantities $1 and 
"^•2 (|73|) , ([7^)) in inverse powers of the hard scale Q and then perform azimuthal angular 
averaging over (p^. 

Expanding $2 in twists in the main logarithmic region 

/i^A^p < < q', (89) 

one observes that twist-1 terms cancel, so one has to retain twist-2 and twist-3 terms pro- 
portional K^: 

$2 = — - 90 
The analogous decomposition for $1 reads 

•^-i = ^4 ^ =r^ 91 

Note that the last term does not contain k^. However, one must track it because it will be 
important in double helicity flip amplitudes. 



4.2 Twist expansion for S wave type mesons 



With the aid of this decomposition one obtains: 
for amphtude L ^ L 



2m 



for amplitude T ^ T 



Q Q 



+ 



2m 



+ 



1 - 2z) 



M 



M + 2m 



M + 2m 2 
2Je k 



Az{l - z) M + 2m 



-k'ieV*)'^ 



Q 



A:le-y*)^--=j(eA)(y*A) 1 



6^2(1 - 2zf 



(92) 



(93) 
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This amplitude is naturally split into s-channel helicity conserving and double helicity flip 
parts 



m' + 2kHz^ + il-z)'') + 



m 



2Q 



M + 2m 



1 - 2zf 2m 



1 + 



4z(l - z) M + 2m 
(94) 



Finally, single spin flip amplitudes are 



It-.l = -2Mz{l - z){l - 2zy^{eA) 

Q 



1 + 



(1 - 2zY 2m 
4:z{l- z) M + 2m 



IUt = -2Qz{l-z){l-2zY^{V*A) 



2/?2 2A;2 M 



Q— 'Q^M + 2m 
4.3 Twist expansion for D— type vector mesons 



(95) 



Here we will need to track higher-twist terms. It will turn out later that leading contributions 
vanish, so twist-3 terms will be crucial for our results. 
For amphtude L — > L one has 



IE-.L - -QMz{i - z) - 

For T —>■ T amplitude, one obtains 



Am 2\ 



4If 



(96) 



9 2/£ 




2te -2 


m ^ 


i'-w] 




Q 




Q 



+ 



2k' 



-4^(1 - z)p' + (1 - 2zym{M + m) < —k\eV*) - -^{eA){V*A) 



Q 



2g 



QKHl-2zf' 

1 + ^^2 

Q . 



2ie 



4P 



-m(M + m)-V{eV*)^ U " =T 



Q 



Q 



(97) 



Note that we kept track of all terms oc Again, one can separate out s-channel helicity 
conserving and double helicity flip parts: 



7|Lt(A, = Av) = (e-y*)=j 



2p2 [ + 2^2 _ 4^2^\ _ + ^)^2 / ^ _ 4^ 



Q 



Q 



tD ^^ XX fr2 4m 2\ /-AN/T-T* AN 6/?2(l-2z) 
/|LT(A, = -Ay) = (k'-—plj{eA){V*A)-=^i^l + 



(98) 
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Finally, single helicity flipping amplitudes are 



-M(l - 2z) 



-2(1 - 2;z)(eA)K2 /r9 4m ,\ , dn'^k^ , 



K 



—,{l-2zyMieA) 



r2 4m 2 \2P 
A; — m M + m — rj- 



-SQz{l -z){l- 2z)2(p2 + + mM)^P(F*A) 



(99) 



4.4 Final results for S wave mesons 



In order to grasp the major features of various S and D wave amplitudes, further simpliflcations 
can be achieved if one neglects spherically non-symmetric arguments of and gluon density. 



First we rewrite general expression (|82D in the more convenient form 



A{x, Q , A) = -ts — / d P77^(p ) / ^asJ" ■ I[i V) 



27r^ 



M 



(100) 



In this expression everything except for integrands 1(7* — > V) is spherically symmetric, thus 
making it possible to perform angular averaging over fip in these integrands. 



4.4.1 S wave: Vt^ averaging 



Here all the calculations are rather straightforward. In the non-relativistic case one can ev- 
erywhere put z — > 1/2 ; M = 2m = my. The resultant integrands are: 



I%T ^T-\^ = Xv) 

I^{T L) 
I^{L T) 



8QM 



(Q2 + M2)2 

8M2 ^2 
(Q2 + M2)2'' 

16 (eA)(V^*A) 
y (Q2 + M2)2 



1 + 



96 



K2p2 



5 M2(Q2 + M2) 



64 M(eA) p2 

T (g2 + M2)2M2 

512 Q{V*A) 



K 



15 (g2 + M2)2 M2(g2 + M2) 



(101) 



Note several things: since the accurate IS* wave differs from the naive 7^ spinorial structure 
only by relativistic corrections, one would obtain the same results in the case of naive qqV 
vertex. The only difference would be only extra factor 2 for the L ^ T amplitude (which is 
higher-twist amplitude, anyway). 
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Thus the only thing left is |p| integration. Note that these amplitudes are naturally 
expressed in terms of decay constants. Indeed, in the extremely non-relativistic case expression 
lOD turns into 

/^^) = ^/rf^P^5 ^ jd'v^s = ^f'' (102) 



27r3 



4.4.2 S wave: answers for L ^ L up to differential cross section 

Here we would like to digress and for the sake of logical completeness show how one obtains 
the final result for the differential cross section with the example of L — L amplitude. If 
needed, the same can be done for the other amphtudes, so we will do it just once. 
One has: 

3 4 -8Qmv fQ^ dn^ dG{x,K^) 1 
^^""^ Gix,Ql)a,(Ql) exp(-i53ipA2) J d'p i^s 



{Q 



exp(--53ipA2)G'(x, Qo)as(Qo)^/i 



m 



2 ^2 

v) 



V 



With this result, (|100|) becomes 

/ycv A/47raem 647r^ Q G ■ as 



A 



is- 



27r2 



3 niy (Q^ + 



327r2 Q 

IS— CvJV ■ y^T^aen 



■G ■ asexp(-|S3ipA' 
3 mv (g^+^P 

The expression for the differential cross section reads it = A^): 

UP 



da 
~dt 



647r3 g2 



9 my 



(cvfv) 47ra 

en 



(g2 + m'^M 



Finally, one can express this cross section through r{V 



V) 

see 



da 
H 



G^a;^exp(-i?3ipt) 



(103) 



4.4.3 5" wave: the other amplitudes 



T — > T, — Ay. 



In the non-relativistic case, this amplitude is readily obtained from the above formulas after 



Q — > my replacement in the numerator of the amplitude (see ( |101| )). This means in particular 
that in this limit 



(104) 
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T — > T, = —Xy. 

This amplitude is very interesting because of the competition of two very different terms — 
soft and hard scale contribution. Indeed, integration over gluon loop gives 

G(a:,^^) 96 G{x,q')p^ 

We see that the soft contribution turns out to be of leading twist, while the pQCD contribution 
is of higher twist. This observation was first made in for the naive type of qqV vertex; here 
we see that it also holds for accurate 5* and D wave vector mesons. 

Since at moderate both terms can be comparable, we cannot neglect soft scale con- 
tribution, which means that we failed to compute the double helicity flip amplitude within 
pQCD approach. 

T -y L and L ^T. 

In the case of heavy quarkonia these single spin flip amplitudes have suppressing non-relativistic 
factors. Besides, the amplitude L — > T is of twist 3, which is another source of suppression. 
Their ratios to A{T ^ T) = A{T ^ T; = A^) read 

A{T ^ L) 8(eA) A{L ^ T) 64 QiV*A) 

The model dependent quantities W2 and W4 are defined via 

1 Jd^pp^i/js 1 Jd^pp^ips 

^■i = — — r:^ — 1 — ; ^4 = — — r-r^ — 1 — • (i07) 

Within the oscillator ansatz (Hsl) their values are 



4.5 Final results for D wave 

This case is much more tricky. It turns out that the leading terms in integrands pro- 
portional to m^lpp cancel out after angular averaging, so that many new terms, including 
higher twist terms come into play. This cancellation is in fact quite understandable. Indeed, 
in the very beginning we showed that vertex 'u'7^m contains both S and D waves, with D wave 
probability being suppressed for heavy quarks due to non-relativistic motion. This means in 
particular that the photon couples to qq pairs sitting either in or D wave state. However, 
at the other end of the quark loop we have a vector meson in pure D wave state. Therefore, 
the largest items in (75+d|---|^d) cancel out due to S-D orthogonality. 

4.5.1 D wave: Vlp averaging for L ^ L amplitude 

If we limited ourselves only to the leading p^/m^ terms, we would get 



dn^ [k'-—pij=A7r. (^-p^ - 2 . -p^ ) = 



which is the manifestation of S-D orthogonahty. Thus, we see that p^/m^ terms vanish after 
angular averaging. Therefore, one has to be extremely careful now and must take into account 
all possible sources of p"^/ m*^ terms. To do so, one has to perform the following averaging 



(4.(1 - .) . . (p - l^p; 

The detailed evaluation given in Appendix y results in 

1 4p^ / M2 




(109) 



M2 



where Qq = m? + One can now again express the integral over quark loop through the 

decay constant (see (^iD): 

(110) 
(111) 



to give 



(g2 + ml ^2 



I' 



15 

Comparison with L L amplitude reveals that 



/ 




(112) 



4.5.2 D wave: the other amplitudes 

For the helicity conserving amplitude one has to repeat the same averaging procedure, which 
is again fully elaborated in Appendix y. The result can be written as 




(113) 

In the case of double helicity flip amplitude we again have contributions from soft and hard 
scales with the same hierarchy of twists. Namely, 

G(x,/x2) 96 G(x,g')p2 

/x2 7M2(g2 + M2) ^''^^ 

so we again conclude that for the case of double helicity flip amplitude soft contribution 
invalidates our pQCD analysis. 

In the case of single spin flip amplitudes, no dangerous calcellations among leading terms 
arise. Before giving a list of amplitudes, we wish to emphasize that in the case of D wave 
mesons there is no non-relativistic suppressing factors like W2 and tU4 defined in (|107|) . This 
means that for moderate momentum transfers helicity non- conserving amplitudes are abso- 
lutely important for that case of D wave mesons. 
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4.6 S wave vs. D wave comparison 

Wc would like to present our final results in the form which stresses the remarkable differences 
between S wave and D wave amplitudes. Below we give a table of the ratios 



Pij 



AS(i ^ j) f(D) 



(115) 



for helicity conserving and single spin flipping amplitudes. Double spin flip amplitudes are 
not given due to the presence of incalculable non-perturbative contributions. 




(116) 



Thus, we note several things. First, the abnormally large higher twist contributions to D 
wave amplitudes are seen here as terms oc my/{Q'^ + my). They even force the opposite sign 
of L — > L amplitude in the moderate domain. Second, we see highly non-trivial and even 
non-monotonous dependence of {All/ AttY ratio, which will lead to the presence of a 
dip in experimentally measured ai j oj- for D wave meson production. Finally, we must stress 
that in the case of D wave mesons there is no non-relativistic suppression for single spin flip 
amplitudes as it was in S wave mesons. This leads us to a conclusion that s-channel helicity 
is strongly violated in the case of D wave meson production. 



36 



Spin density matrix for dijfractive p at HERA 
(Schilling- Wolf conventions) 
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Figure 10: The p meson spin density matrix: comparison of recent HERA data with our 
predictions for IS* state. 



5 Numerical analysis 

In this section we will present some results of the numerical investigation. These results bear 
mostly illustrative purpose. Undoubtedly, they do not exhaust all interesting issues and leave 
much room for further numerical analysis. 

Everywhere in our analysis we used oscillator type wave functions (|45|) with Rd = Ris- 
There are several parameters to vary in order to achive the most correct description of the 
vector mesons. Aside from the very ansatz for the wave function, which is now fixed, we 
also have quark mass m and radius R as free parameters. We stress that these are purely 
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(A(LL)/A(TT) )**2for rho system 
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Figure 11: Ratio [A{L L)/A{T 
a function of Q^. 



T)Y for S and D states in p system at fixed x = 0.01 as 



phenomenological parameters; in particular, it would be a mistake to assign values of several 
MeV to the quark masses. In principle, these quark masses can even vary from the ground to 
excited states for the same quarkonium. 

We used normalization conditions (|32D and ( plD to determine normalization constants and 
chose a specific pair m, R to have the decay constant coinciding with the experimental values. 
To avoid unnecessarily many free parameters, we used the same values for 5* and D wave states. 
Specifically, for the p meson we took m = 0.3 GeV and R = 4.0 GeV~^ (which corresponds 
0.8 fm). For charmonium the value were m = 1.3 GeV and R = 1.5 GeV~^ =0.3 fm. 

Because all quantities displayed in figures are ratios of one amplitude to another, these 
results are not sensitive to the particular gluon distribution parameterization. Since we leave 
double helicity flip amplitudes aside, we can now switch to the integrated gluon density 



So, we used the integrated form of GRV parameterization with a linear extrapolation of gluon 



density in the region Q < \ GeV so that G{x, Q 



when Q 



0. 



In fact, in an 

accurate numerical analysis low region should be treated more carefully. Indeed, parton 
distributions do not evolve in this region, however, they can be parameterized with aid of low 
energy experimental data. 

In Fig.^ we compare our calculations for p vector meson density matrix with recent 
experimental data from HERA [^. We see that our calculations describe data reasonable 
well, including parameters connected with s-channel helicity violation. 



In Fig.0 we show ratio {All I AttY (which would be olIot if the s-channel helicity were 
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\A(TL)/A(TT)\ 



rho system 

0.4 ^ ' 




0.05 



10 100 1000 

Figure 12: The absolute value of ratio A{T L)/A{T T) revealing the relative magnitude 
of T ^ L spin flipping amplitude at fixed x = 0.01 as a function of Q^. The sign of this ratio 
is — for S wave and + for D wave states. 

conserved) for both S and D wave states of p and J/ ip system. A striking difference is observed 
both at small and large Q^. 



In Figs.|l2|, [T^ we present the relative magnitude of helicity non-conserving amplitudes 
again for S and D states in p and J/ip system. It is interesting to note that T ^ L ampli- 
tudes have opposite signs for 5* and D wave mesons. We underline the large helicity violating 
amplitudes in D wave states, both in non-relativistic J / ip and fully relativistic p meson. The 
size of the effect reaches 35% for p system and 13% for charmonium. Thus, it becomes now 
quite obvious that in the case of D wave mesons no helicity conservation holds even at the 
qualitative level, in accordance with our previous conclusion from the heavy quarkonia analysis. 

Finally, we analyzed the possibility of S — D mixing in the physical p meson. Namely, we 
introduced mixing angle 9 as 

Pphys = cos6'|15') + sm9\D) (117) 
and plotted ratio (All/Att)'^ in Fig. 14 for mixing angles 6 = ±20°. The results are strikingly 
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-A(LT)/A(TT) 

rho system 

0.3 ■ r 




^ • — 

10 100 1000 



2**2, GeV**2 

Figure 13: The ratio —A[L — > T)/A{T — > T) showing the relative magnitude of L — > T spin 
flipping amphtude at flxed x — 0.01 as a function of Q^. It is seen that this transition is of 
higher twist. 

sensitive to the presence of D wave admixture in the physical p meson; further investigation 
of this phenomenon is needed. 
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Figure 14: The impact of possible S — D mixing in the physical p meson on the ratio [A{L — > 
L)/A{T T)]2 at fixed x = 0.01 as a function of Q^. 

6 Conclusion 

In this work we brought under an accurate analysis the issue of internal vector meson spin 
structure in its diffractive production in DIS. In fact, we developed an essentially new way to 
study the internal spin structure of a vector meson. 

To do so, we elaborated a formalism of treating pure S and D wave vector mesons in such 
processes. For the first time, the full sets of helicity amplitudes in both cases were calculated. 
The exact results (|87D , (^) were derived within pQCD leading log approximation and 
then analyzed analytically in the limiting case of heavy vector mesons ( |103| ), (|104|) , ( 106 ), 
(|116|) and numerically in the case of p system and charmonium (Figs. 10-14). 



In particular, for the s-channel helicity conserving amplitudes we observed: 

• remarkably different Q^-dependence of S and D wave type amplitudes, providing thus 
a way to discern S and D wave states that are indistinguishable at e~^e~ colliders; 

• dramatic role of higher twist contributions to the D wave vector mesons, which even 
forced sign change for L — > L amplitude and led to non- monotonous (Jl/ctt ratio; 

• in the case of light vector mesons we found that D wave amplitudes are absolutely 
comparable with S wave, making, for example, the physical p meson very sensitive to 
the possible presence of D wave admixture. 

Our results for hehcity violating amplitudes include: 

• observation that they are not negligible for light vector mesons of S wave type, but 
become suppressed for heavy S wave mesons; 
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• very large helicity violating effects for D wave vector states, which do not get suppressed 
even in the case of heavy quarkonia; 

• opposite signs for the largest spin flip amplitude T L for S and D wave vector mesons; 

• confirmation of the soft dominance of the double spin flip amplitude in the case of 
accurate S and D wave states. 

Though we performed some numerical analysis, there is much room for further investi- 
gation. In particular, more extensive analysis is needed to check the model dependence of 
the numerical results obtained. It is also desirable to give definite numerical predictions for 
diffractive DIS experiments. 
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A Denominator evaluation: details 



Below we give a detailed derivation of Eg. (pBD . The major guideline will be again analysis of 
pole positions and setting some of the propagators on mass shell by taking appropriate residues. 

We first start with the case of strictly forward scattering. (Though it corresponds, of 
course, to forward elastic 7* — >■ 7*, not '-f* —>■ V scattering, it will not have any effect on the 
expressions we derive here.) So, the integral to be calculated (see Fig.^]) is 

^ ^ [(/c — k)^ — rn? + it] [/c^ — m? + it] [{k — qY — m? + it] [{k — q — kY — rn? + it] 

(118) 



[k^ — ji^ + it^lip — kY — + it] 
With Sudakov's decomposition (|47|), one can rewrite all 7 propagators as: 
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;ii9) 



After an accurate analysis of position of all poles on complex a, /?, y planes, one finds that the 
only non-zero contribution survives if one uses line (4) to extract a, line (2) to extract y and 
line (7) to extract 13: 

(k — kY + m? lY + m?' ^ + iM^ , , 

a = -y + x+ ^ ' ; y = ; (3 = -. 120 

(1 — z)s zs s 

Note that resultant a,l3,y ^ 1. The remaining z-integration goes from to 1. Therefore, 
integral ( |118| ) turns into 

27ri\3 dz 1 1 

' W 'io ~T^' (5)(6)(1)(3) ' ^ ' 

where the numbers in brackets denote the corresponding propagators. With aid of (|120|) , these 
propagators become: 

(5,6) -aPs - (k^ + fY) ^ -{k^ + fY) 

(3) y{l - z)s - xs{l - z) - {JY + m^) = ~-jY + rrY + z{l - z)Q'^] 

(1) -{y + a- x)zs - [{k - kY + = " '^)^ + + z{l - z)g^](.122) 



Here we used a <^ 1. Substituting these expressions into ( |121| ), one gets 

[]' . f _± \ n23) 

) Jo {k^ + fi^Y [k^ + + z{l ~ z)Q^][{k - kY + + z{l - z)Q^] 



27ri\3 
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Note that we replaced the lower integration limit by zero, knowing that there is no large end- 
point contribution. 

In the above analysis we took momentum transfer A = for simplicity. Now let us see 
what non-zero momentum transfer will change (see Fig.@). 

First of all, now A will appear in gluon propagators. However, condition ( pBD allows us 
to neglect it within diffraction cone. Then, A will enter quark loop propagators. Here the 
presence of momentum transfer (especially, the longitudinal one) is very important. Indeed, 
within our notation of momenta, the propagator pair (1,4) will yield the same contribution 
as before 

kl + m^ + z{l - z)Q^ , 

where ki is transverse momentum flowing through the photon vertex. On the other hand, 
propagators (2, 3) will give now 

fc^ + — z{l — z)'my , 

which together with vertex factor F gives rise to the vector meson LCWF. Combining every- 
thing together, one now has 

r F 

/ dy dz da d(3 -r—rz r 

J [all propagators] 

_ / 2^1^^ r dz p 1 1 
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B Helicity amplitude technique evaluation 



Here we give the derivation of expression for traces of the following type 
Sp |e (A;4 + m) p' {k^ + m) V* {k2 + m) p' {ki + m)| 



(125) 



in full detail. Though one can calculate this trace covariantly, a particularly convenient way 
to do so is given by light cone helicity amplitude technique [0. We emphasize that both 
ways are absolutely equivalent. In the helicity amplitude approach, we recall that all fermion 
lines in ( |125| ) can be taken on mass shell (see detailed derivation of LCWF normalization in 
Sect. ^75|) and decomposed into a sum of light cone helicities 

(k + m) ^ UxUx for quark lines; 

A=± 

{k + m) = — [{—k) — "W-] ^ — ^ VxVx for antiquark lines. 

A=± 

Since the specific choice of this decomposition does not affect the final result, we are free to 
take the most convenient choice of spinors (see for details), namely, 

1 



u{p, A) 
v{p,X) 



V2p^ 
1 



where 



(V2p'^ + l3m + dp 
\/2p^ — (3m + dp 




1 

V J 



^^^^ = 71 



x(T) 


A = 


+1 




A = 


-1 


xii) 


A = 


+1 


x(T) 


A = 


-1 


( ° ^ 






1 













V -1 ) 







(126) 



It is a straightforward exercise to show that they obey usual normalization and completeness 
requirements: 



u{p, \)u{p, A') = -v{p, \)v{p, A') = 2m5xx' ; 

u{p, X)u{p, X) = p + m; ^ v{p, X)v{p, A) = p — m . 



Thus, expression ( |125|) turns into 
2s2./W = 



— 1)^ ^ uev ■ vp'v ■ vV*u ■ up'u ■ 



(127) 



(128) 



Though this expression is written for Diagr.|^, one can immediately redo the same transforma- 
tion for the other diagrams. So, in the following we first list the light cone helicity amplitudes 
to be used later. Then, we derive the amplitudes for photon and vector meson vertices uev 
and vV*u respectively with appropriate polarization vectors. Finally, we give a simple rule of 
how to include 2 gluon legs attachments (i.e. up'u and vp'v) and write the final expressions. 
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B.l List of helicity amplitudes 



Here we list various amplitudes of u')^^{p)...u\^{q) type various gamma matrix structures in- 
serted between light cone spinors. This notation should not be confused with our previous 
usage of p and g; here and only here these momenta are assigned to the fermion lines as shown 



at Fig. 15 



For a derivation of these amplitudes see [|T7 . 



a) 



q 

P 




c) 



b) 




d) 



Figure 15: Four types of transitions for which the light cone helicity amplitudes are given 
below. The dot with label A indicates the spinor, whose helicity is used as A in the table of 
amplitudes. 

Notation: 

as usual, each vector a should be understood as a±, 

momenta p and q are constituent momenta shown at relevant figures; 

For A = ±1 we define 

a(A) = —Xa^ — iay . (129) 

'Cross product' is defined as 

[ab] = ajjy - bxtty . (130) 

(This is not a true cross product of two vectors, this is just a handy notation.) Moreover 
every matrix element should be multiplied by common factor ^p^q^. 



B.1.1 Straight line elements 



Fig. 15a: 



ux{p)...ux{q) 



u.x{p)-ux{q) 



7 
1 



+ 



^ {m? +pq + i\\pq\) 



P+Q+ 



" ' p+ ^ g+ ^ \ p+ q 



m 



(p(A) - g(A)) 

—■ma{X) ( 

\P+ q+j 



(131) 
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Fig.||b 

a ■ 7 

B.1.2 Vertex lines 



vx{p)...vx{q) 



vx{p)-.v-x{q) 




{w? + pq — i\ [pq] ) 



m 



ap , dq , „- x f[ap\ [aq] 

p+ ^ q+ ^ \ p+ q+ 



p+q+ 



ma{X) 



{p{\) - q{\)) 



1+. 



Fie;. 15c 



vx{p).-ux{q) v^x{p)--ux{q) 



7 

7" 



p+q+ 





(p(A) + g(A)) 



p+q+ 



-m 



+ pq + iA[pg]) 



a ■ 7 ma 



+ - + - 



iX 



p+ q+ 



p+ q-\ 



Fig.lSd 



7 

7' 



ux{p).-vx{q) 




ux{p)-.v^x{,q) 
2 



m 



(p(-A) + g(-A)) — ^ — i-m^ + pq + i\[pq\) 



a -7 -"^«(-A) + ^) 







B.2 Photon vertex amplitudes 

Notation is given in Fig.^ We start with transverse photon. 

Equal qq helicities give 



,Jz{l-z)q+ 



^(1 - z)q+ 



-m)e(— A) 



me(— A) 



Opposite qq helicities give 



^(1 - z)q+ 

^(1 - z)q+ 
1 



(efci)(l - - (efci)z - iA {[eki]{l - z) + [eki]z) 



{eki){l-2z) -iX[eki] 
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In the case of scalar photon 



u'eov = ti'^(g+7_ + xp_'j+)v 
the same hehcities give exactly zero while opposite helicities result in 



zil - z)q+ 1 



z{l - z)q+ Q 
1 1 

z(l - z)Q 



— i—rn? — A;?) + x2z(l — z)p_q> 



m^ + kf- z{l - z)Q' 



(135) 



B.3 Vector meson vertex amplitudes 

This case is more tricky due to the nonzero vector meson transverse momentum A. We start 
with the transverse vector meson polarization: 

^rr* 2(\7*A) 
V* = -^V* + ^ '-P-1+ ■ 



The same qq helicities give again 



_ ^z{l-z)q+ 
^(1 - z)q+ 

while opposite helicities give 

\jz{l - z)q+ 



mV*{\) 



z{l-z) 



mV*{\) 



z(l — z)q 



[{V*h)z + {V*k2){l -z)-i\ [[V*k^]z - [V*k2]{l - z 



p^2z{l — z)q^ 



{{V*k){l - 2z) + iX[Vk]} 



'z{l-z) 

Here we used definitions and properties (see also Fig.||): 
k2 = k + zA; k3 = -k+ {1- z)A 

{1 - z)k2 - zh = k ; {1 - z)k2 + zks = {1 -2z)k + 2z{l - z)A : 
ik^k^) = -p + (1 - 2z){kA) + z{l - z)A^ 

2 

+ A\ 



9 1*9 kl + m? ki + m? kP' + m? 
M + A = — — ' 



z ^ [l-z) z[l-z) 
For the longitudinal vector mesons one has for equal quark-antiquark helicities 



^(1 - z)q+ 1 



z(l - z)q+ M 



-mA(A)+m^[A;2(A) + A;3(A)] 







(136) 



(137) 



(138) 



(139) 
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and for opposite helicities 



2^(1 - z)q+ 1 
z{l - z)q+ M' 
A2-M2 



-[{Kh)z + (A4)(1 -z)+ + ^[-m^ + (^fcg) + zA[44] 

9+ 



+ 



2z(l - 2)M 



(140) 



B.4 Final trace calculation 



Once we have the expressions for vertex amphtudes, the rest is done quickly. We first note 
that each gluon vertex attached to the lower or upper line gives factor 



2zqj^ ■ = sz ; 2(1 — z)q+ ■ p_ = s{l — z) 

correspondingly. 

So. let's start with T —* T amplitude and calculate it for Diagr.(c) at Fig.^. 
Equal qq hehcities give 



'1411 



s(l — z) ■ sz 



Summing over A gives 



zil-z) 



--mV*{X) 



=me(-A) = -sWe(-A)\/*(A) 



z(l-z) 



The opposite hehcities yield 



[{V*k){l - 2z) + iX[V*k]\ [{eki){l - 2z) - iX[eki] 
Summing over helicities and making use of identity 

[a6][c(i] = {ac){bd) — {ad){bc) 

one obtains 



2s' 



{V*k){eki){l - 2zf + {eV*){kki) - {ek){V*ki) 



Since we factored out 2s^ when deriving (|6^) , we finally get 



{eV*){m^ + kki) + {V*k){eki){l - 2zf - {ek){V*ki 



(142) 

(143) 
(144) 

(145) 
(146) 

(147) 



We included factor (—1) since in this diagram we have one antiquark propagator. 

An important observation here is that all other integrands, namely J^"-* -{l — z)/ z^ I^^^ , I^'^^ ■ 
z/{l — z) give absolutely the same result (with their own definitions of ki of course). The only 
thing one should not forget is that diagrams (a,d) enter with sign '— ' while diagrams (b,c) 
enter with sign '+': 



:jH ^ j(fe) ^ j(c) 



l~z 
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For L L amplitude one has immediately 



41l = -^h' + z{l - z)Q'] . ^2z{l - z)M' (148) 

In fact, using simple relation 

m' + kf- z{l - z)Q^ ^ ^ ^ -2z{l - z)Q^ 
-m? + kl + z{l- z)Q'^ m'^ + kl + z{l- z)Q'^ 

and noting that all unity terms will eventually cancel out in (^), one can rewrite (|148|) as 

IlIl = -^QMz^{l - zf . (150) 

For T —>■ L amplitude one has 

I^l^ = 2z{l- z)M{eh){l-'2z) (151) 
and for L —>■ T amplitude one has 

= l[m2 + kl- zil - z)g2](i _ 2z){V*k) . (152) 

The same transformation as in L — > L amplitude, leads to 

I^^l^ = -2z{l - z)Q\l - 2z){V*k) . (153) 
Note that in the last three amplitudes only opposite qq helicities contributed. 
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C Averaging over Qp for D wave states 

Before doing this, let's have a table of useful relations: 



{k'P) = ^^p' m) = ^^p' {pIpI) = ^^p' (154) 



Finally, remember that = |(1 — 2z)^M^. 
One has to perform the following avegaring 



2 



(4.(l-.).i,.(P-i;^ri).(i_«|> ,55, 

Note that all factors should be carefully examined; all four do contribute to the final answer. 
Decomposing Q as 

Q' = m' + zil - z)Q' = ^2 _^ 1q2 _ 1 (1 _ 2zfQ' ^Q]- ^Q' , (156) 

one gets 



_4 - _4 .1 + 2^^^ . (157) 



So, omitting Qq^, one has 



APJ \ hpQ'J V A/-V g; 
With the aid of ( |100| ), one obtains 

2 2_4m 2\_A_ 4/2 6 A I o 4/l__6_A_A ^ ('A _ 

s"^ sm"^ y U5 157 M2go'^ u5 15; qI^ u5 15; 

2 4p2 16p^ Sg^p^ 16p^ 
3*^M + 2m 15M2 ISM^Q^ 15Qo 
4p4 ^ 4 8 g2 16 M2 \ 



3M2 \^ 5 5Q2 + M2 5Q2 + ]Vf2^ 



15M2 \^ Q'^ + M'^J 

For the helicity conserving T —>■ T amplitude one has to repeat the same averaging proce- 
dure. The expression to be calculated first is 



1 + 2 



2p2 + 2A;2 - Ak'^^ - m{M + m)k^ 1 - =^ - 2A;M A;^ - — 



g / V Q 

52 



^ 3 ^ M^qI\2 3 157 3^ 15^ Qo 

4_A 4,8 4 16 4,8 4 

3^ 15^ 3^ 15M2 + Q2P + i5M2 + g2P 
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